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Abstract. We prove Lieb-Thirring inequalities for Schrodinger operators with a 
homogeneous magnetic field in two and three space dimensions. The inequalities 
bound sums of eigenvalues by a semi-classical approximation which depends on the 
strength of the magnetic field, and hence quantifies the diamagnetic behavior of the 
system. For a harmonic oscillator in a homogenous magnetic field, we obtain the 
sharp constants in the inequalities. 



1. Introduction and main result 

Lieb-Thirring inequalities |LiTh] provide bounds on the sum of negative eigenval- 
ues of Schrodinger operators in terms of a phase space integral. In this paper, we 
are interested in two-dimensional Schrodinger operators Hb + V with a homogenous 
magnetic field of strength B > 0. Here 

. d Bx 2 \ 2 ( . d Bx 1 
H B = -i-z h — — + \ -i 



dx\ 2 J \ dx 2 2 

is the Landau Hamiltonian in L 2 (M 2 ) and V is a real- valued function. The Lieb- 
Thirring inequality states that 

Tr(H B + V)_<r 2 {2n)- 2 [[ (\p\ 2 + V (x)) _ dx dp (1) 

with the (currently best, but presumably non-optimal) constant r 2 = 7t/a/3 from 
[DoLaLo] . Physically, the left side is (minus) the energy of a system of non-interacting 
fermions in an external potential V and an external, homogeneous magnetic field of 
strength B, whereas the right side is — r 2 times a semi-classical approximation to that 
energy. 

Physically, one expects the system to show a diamagnetic behavior, that is, to have 
a higher energy in the presence of a magnetic field. This is however not reflected in (JT]), 
which has a right hand side independent of B. We refer to [Ft] for further references 
and a survey over this problem. Our goal in this paper is to obtain a bound similar to 
(CD), but with a more refined semi-classical approximation which takes B into account. 
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The approximation we propose is 



(2) 



B f 

— V/ {{2m+l)B + V{x))_ dx. 

Z7T z — 4 /ma 

This quantity reflects the diamagnetic behavior since 

I ((2rn + l)B + V(x))_dx<-±- If (\p\ 2 + V (x)) ^ dx dp (3) 



for every V^. Inequality (J3]) follows (even before the ^-integration) from an easy convex- 
ity inequality (see Lemma [121 below) . We also note that when B — > 0, by a Riemann 
sum argument, the quantity ([2]) approaches 

(An)- 1 [ dE [ (E + V(x))_ dx = (2tt)- 2 [ [ (\p\ 2 + V(x)) __ dx dp , (4) 

JO Jr 2 J JR 2 xR 2 

which is the 'usual' phase space integral. 

While the right side of fll]) (up to the constant r 2 ) has the correct limiting behavior 
when a small parameter h is introduced, it is not useful in the coupled limit B — > oo 
and h — > 0. This limit is physically relevant, for instance, in the study of neutron stars 
|LiSoYn] . The magnetic quantity (T5]) reproduces the correct behavior in this regime. It 
is remarkable that this asymptotic profile is, indeed, a uniform, non-asymptotic bound. 
This is implicitly contained in |LiSoYn2] who use, however, only an approximation of 
02]). Our first result is 

Theorem 1. For any B > and any V on M 2 one has 

b 00 r 

Tr{H B + V)_ <p 2 — V / ((2m + l)B + V(x))_ dx (5) 

with p 2 = 3. 

Hence, up to the moderate increase from r 2 = 7t/a/3 ~ 1.81 to p 2 = 3, we have 
found a magnetic analogue of ([I]) which reflects the desired diamagnetic behavior (|3]). 
An important ingredient in our proof is a method developed recently by Rumin |Ru] 
to derive kinetic energy inequalities; see Subsection 12.11 



Similarly as in the non-magnetic case, one might ask for the optimal value of the 
constant p 2 - By the semi-classical result mentioned above one necessarily has p 2 > 1. 
A first result in this direction was obtained in [FrLoWe] (extending previous work of 
|ErLo Voj ) . where it was shown that if one takes V to be constant on a set of finite 
measure and plus infinity otherwise, then ()5]) holds with p 2 = 1. Our second main 
result is an analogous optimal bound for a harmonic oscillator. 

Theorem 2. For any B>0,Ui>0,U2>0 and /i > 0, inequality ([5|) holds with 
p 2 = 1 for V(x) = bj\x\ + <J\x\ — fi. 
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In particular, letting B — > and using the limit in (J4j) we recover the known bounds 
in the non-magnetic case from |dB} ILa2] . Even though the eigenvalues of a harmonic 
oscillator in a homogeneous magnetic field are explicitly known (Lemma ITUl). the proof 
of Theorem |2] relies on a delicate property of a subclass of convex functions (Lemma IT4l) 
which, we feel, could be useful even beyond the context of this paper. 

Moments of eigenvalues. Using some by now standard techniques we derive a few 
consequences of Theorems [U and El First, following Aizenman and Lieb |AiLi] one can 
replace V by V — p in (|5]) and integrate with respect to p to obtain that for any 7 > 1 

Tr (H B + V)l<p 2 —Y / ((2m + 1)B + V(x))l dx , (6) 
27r ^o^ 2 

where p 2 = 3 for general V and p 2 = 1 for V(x) = <J\x\ + (jj\x\ — P- The restriction 
7 > 1 is necessary, since one easily checks that for < 7 < 1 there is no constant p 2 
such that ([6]) holds for all potentials V. Restricting ourselves to the quadratic case we 
shall show in Subsection 13.41 

Proposition 3. For any < 7 < 1 there are B > 0, p > and uji = u 2 such that for 
V(x) = uj\x\ + uj\x\ — p one has 

B 00 r 

Tr(H B + Vy_>— V/ ((2m + 1)B + V(x))l dx . (7) 

2VT L — ' /ra>2 

Our counterexample appears in the limit ujj/B — > (with p/B = 3 fixed). 

Three dimensions. Next, we shall show that our bounds for d = 2 can be applied to 
deduce analogous bounds for d = 3. This argument is in the spirit of the lifting argu- 
ment from |Lal4 \La2\ ILaWej . We denote by H B = H B — the Landau Hamiltonian 
in L 2 (R 3 ). 

Corollary 4. For any B > and any V on R 3 , one has 

B 00 f 

Tr(H B + Vy<p 3 —-J2 /. ((2m + l)B+pl + V(x))_dxdp 3 (8) 
with p 3 = a/3 71"- 

Proof. From the operator-valued Lieb-Thirring inequality of [Do LaLo] we know that 

Tr {h b + V) _ < ^ jj^ Tr LHm (H B + pg + V(; x 3 )) _ . 

Inequality OH]) is therefore a consequence of Theorem [TJ □ 

For the harmonic oscillator we have 

Corollary 5. For any B>0,ui>0,u 2 >0,u 3 >0 and p > 0, inequality (jSJ) holds 
with ps — 1 for V(x) = uj\x\ + uj\x\ + uj^x^ — p. 
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Proof. We denote by Ej the eigenvalues of the one- dimensional harmonic oscillator 
H = - JSj + u%xj. Then, since H B + V = (H B + V) ® I + I <g> H with V(x 1 ,x 2 ) = 
uj\x\ + io\x\, we have 

Tr i2(K3) (H B + V^ _ = J2 T r L 2 (R 2 ) (#b + V + Ej - p)_ . 

i 

According to Theorem [2] (which trivially holds for p < as well) , this is bounded from 
above by 

<x> 



B . ^ r 

— / / / ((2m + 1)B + V{x\, x 2 ) + — dx% dx2 

Z7T Z — 4 Z — 4 /ra2 

oo „ 

} j I Tr L 2 (R) (H + (2m + l)B + V(xi, x 2 ) - ctei dx 2 



2?r 

m=0 ' 



Next, we shall use that H satisfies a Lieb-Thirring inequality with semi-classical con- 
stant 0H1 ILa2j . that is, for any A G R, 

Tr L 2 (R) (H - A)_ < — j I (p 2 3 + u\xl - A) _ ds 3 dp 3 . 

(This can also be seen from Lemma 1121 and recalling the explicit form of the eigenvalues 
of H.) It follows that for every fixed (xi,x 2 ) 

Tr i2(R) (H + (2m + l)B + V{x u x 2 ) - p)_ 

<7T II {pl + ^lxl + {2m + l)B + V{x 1 ,x 2 )- dx 3 dp 3 , 
which proves the claimed bound. □ 

Remark 6. The previous proof shows that ([S]) with p 3 = 1 is valid for more general 
potentials V(x) = V(x\,x 2 ) + v(x 3 ), where V(xi,x 2 ) = uj\x\ + u 2 x 2 and where v is 
such that Tr L2(R) (-^ + v(x 3 ) - A)_ < ^ jJ RxR G» 3 + v(x 3 ) - A)_ cfe 3 dp 3 for all A. 

A similar argument as in the proofs of Corollaries H] and |5] (based on the operator- 
valued Lieb-Thirring inequalities of |HuLaWe"| ILaWej ) shows that for general V one 
has 

Tr (h b + vY < P3i7 —_ V / ((2m + i) B + p 2 3 + V(x)) 7 dxrfp 3 (9) 

with p 3j7 = 6 if 7 > 1/2, with p 3j7 = 7rv^3 if 7 > 1 and with p 3 . 7 = 3 if 7 > 3/2. 
Moreover, in the special case of V(x) = uj\x\ + uj\x\ + w 3 x 3 — p, (Q holds with p 3 = 1 
for 7 > 1 and with p 3 7 = 2 (7/(7 + l)) 7 for < 7 < 1. The latter follows from the 
fact (FrLoWej that 

(-£ + ■** - a ) 7 _ £ 2 fc?r) 7 s / L (pl + _ A)! dx:> dn ■ 
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2. Proof of Theorem [T] 

2.1. A kinetic energy inequality. We define a piecewise afine function j : [0, oo) — > 

[0,oo) by 

j(p) = f- (P 2 + (2L + l)r) if p = ^{L + r), L G N„, r G [0, 1) . 

Z7T Z7T 

We note that j is continuous, increasing and convex. One has j(p) = Bp if p < B/(2n) 
and j(p) ~ 27rp 2 if p ^> B. The connection between this function and the right side 
of (jSJ) will become clearer in the next subsection. 

Theorem 7. Let < 7 < 1 be a density matrix on L 2 (M. 2 ) with finite kinetic energy. 
Then 

TtHb^ > 3 / j(p y (x)/3) dx , 

</R 2 

where p 7 (x) = 7(2;, x). 

It is easy to see that 3j'(p/3) > (1/3) j(p) for all p > 0, and therefore we also have 

Tr# B 7> (1/3) / j{p 7 {x))dx. 

Proof. The first part of our proof follows the method introduced by Rumin [RuJ. We 
define j R : [0, 00) [0, 00) by 

2 



Mp) = b p + 2bJ2(vp-^) ■ 

k=l v J + 



We note that is different iable and convex, jn(p) = Bp if p < B/{2n) and jn(p) ~ 
2np 2 /3 if p > 5. We shall first show that 

Tr# B7 > / j fl (p 7 (ar)) . (10) 

</R 2 

In the second part of our proof (see Lemma [8]) we show that ji?(p) > 3 j(p/3) for all 
p > 0. 

For the proof of (1101) we write 

Tr H B i= / Tr(P E 7 )dP = / / p E (x)dEdx, (11) 

where P E is the spectral projection of Hb corresponding to the interval [E, 00) and 
where p E (x) = (P E 'jP E )(x,x). It is well-known that 

(1 - P E ){x, x) = — #{m G N : (2m + 1)B < E} . 

The same clever use of the triangle inequality as in |Ruj leads to the pointwise lower 
bound 

P$(x) > l^)-^#{meN : (2m + l)P<P} 
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Inserting this bound in ( TTUl) we obtain 

r I pB °^ r(ZK+ 

TrH Bl > / / Pl (x)dE + V / 

JR 2 \ JO k=1 J(2k-1 



B ^ y(2fc+l)B / , 

^ ( yp 7 (x) - y — ] J r/.r 



j R (p y (x))dx. 

This completes the proof of (fTOl) and also, by Lemma[8]below, the proof of the theorem. 

□ 

Lemma 8. Jr{p) > 3 j(p/3) for all p > 0. 
Proof. We are going to prove that 

j R (3 P )>3j(p). (12) 

Note that this is an equality for p < B/(6tt). Moreover, since the left side of (TT2|) is 
convex and the right side linear for p < B/(2ti), we conclude that (fT2|) holds for all 
P<B/(2n). 

Henceforth we shall assume that p > B/(2ir) and we write 3p = (B/2ir)(K + s) 
with K G N and s G [0, 1). If K = 3L + m with L G N and m G {0, 1, 2}, then the 
lemma says that 

AT + s + 2 ^ (yWTs -Vkj >3(L 2 + \(2L + l)(m + s)) . 
fe=i 

We expand the square on the left side and insert L = (K — m)/3 on the right side. 
This shows that the assertion is equivalent to 

K 

K + s + 2K(K + s) - AVK + s^Vk + K(K + 1) > \K 2 + \Ks + s + R, 

fc=i 

for ifeN and s G [0, 1), where i? = — |m 2 — |ms + m. Since the inequality has to 
be true for any m G {0, 1, 2}, we can replace R by its maximum over these m (with 
fixed s), that is, by (2/3) (1 — s). Thus (fl2"|) is equivalent to 

K 

AK 2 + (3 + 2s) AT - 6VK + s ^ - 1 + s > . 

fc=i 

The proof is straightforward for K — 1 and we may therefore assume that if > 2. By 
the concavity of the square root we have 

Vk + VkTl ^ f k+1 



2 

Summing this from k = l to k = K — lwe get 



< Vidt. 
Jk 



Y^y^k< J Vidt + 



l + VK _ 2K 3 / 2 K 1 ' 2 1 
2 ~ 3 + 6 
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A' 



AK 2 + (3 + 2s)K - 6y/K + s ^ Vk -1 + s 



k=l 



> AK 2 + (3 + 2s)K - y/K(K + s)(4K + 3) + y/K + s + s 

sK((As - 12)K - 9) 

= , h y K + s - 1 + s . 

AK 2 + (3 + 2s)K + ^/AT(A:+-s)(4A: + 3) 

In the quotient on the right side we estimate the numerator from below by — 3sK (AK + 
3) and the denominator from below by AK 2 + 3K + K(AK + 3) = 2K(AK + 3). Thus 
the quotient is bounded from below by — 3s/2, and we conclude that 

K 3R s 

AK 2 + (3 + 2s)K - 6V K + s^Vk - — > VAT + s - 1 - - . 

it=i 

The right side is easily seen to be positive for K > 2 and s G [0, 1), and this concludes 
the proof. □ 

2.2. Proof of Theorem [JJ. In this section we are going to deduce Theorem [T] from 
Theorem [71 We define 

oo 

P(v) :=-— V((2m+l) J B + ^)_ 

Z7T z — ' 

m=0 

for u G R. This is a convex, decreasing and non-positive function. The key observation 
is that this p is the Legendre transform of the function j from the previous subsection, 
that is, 

P( v ) = ™f + V P) ■ ( 13 ) 

p>0 

This can be verified by elementary computations. 

In order to prove Theorem [1] we apply Theorem [7] to get the estimate 

Tr (H B + V) 7 > / (3j(p 7 (a:)/3) + V{x) Pl {x)) dx 

for any < 7 < 1. According to ( {TBI this is bounded from below by 3 L 2 p(V(a;)) dx. 
For 7 equal to the projection corresponding to the negative spectrum of Hb + V we 
obtain the assertion of Theorem [TJ 

Remark 9. Similar arguments show that Theorem [7] can be deduced from Theorem [TJ 
Indeed, since j is convex it is its double Legendre transform. By f|T3|) we obtain 

j(p) = m£(p(v)+vp) . (14) 

By the variational principle and Theorem [TJ we can estimate for any < 7 < 1 and 
any V 

Tri^B7 > -Tr (H B + V)_ + / V(x)p J (x)dx> (3p(V(x)) + V{x)p 1 (x)) dx . 

Jr 2 ' Jr 2 
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According to ( |14l) this is bounded from below by 3 f R2 j(p 1 (x)/3) dx, and this shows 
Theorem [7J 

3. Proof of Theorem [2] 

3.1. The spectrum of Hb + V. The explicit form of the eigenvalues of Hb + oj 2 \x\ 2 
was discoverd in [Fb]. We include an alternative derivation of this result, which is also 
valid in the non-radial case. 

Lemma 10. For any B > andu>i,u> 2 > the operator H b + oj\x\ + io\x\ has discrete 
spectrum and its eigenvalues, including multiplicities, are given by 

B Mt, f) (2* + 1) + a_(f, f ) (21 + 1)) , k, I e N , 



a±K V2) = Ji(l + °l + °l ± ^/(l + °\ + - 4fffa^j . (15) 

Remark 11. It will be important for our analysis below that 

a-(a) a + (a) = a ± a 2 , (16) 

which is easily checked. 

Proof. By means of the gauge transform e - tBx ^ x ^/ 2 W e see that Hb + V is unitarily 
equivalent to the operator 

d 2 f 8 \ 2 

I ■ u — » .22 i , .2 2 



d ^ + ^-i—-Bx 1 j + u 1 x l +u 2 x 2 , 

which, in turn, by a partial Fourier transform with respect to x 2 , is unitarily equivalent 
to 

+ [x 2 — Bxi) z + ~ w 



f) 2 H 2 

° I / D \2 . 2 2 2 ^ 

^ + (x 2 - SxO +w 1 x 1 -a; 2 ^2 . 



2 

t 



After scaling x 2 i— )■ o^a^ this becomes the non-radial harmonic oscillator — A + or Ac 
with the matrix 

'£ 2 + w 2 -5u; 2 N 
-Boj 2 u 2 

The eigenvalues of A are B 2 a + (ui/ B, u 2 /B) 2 and B 2 a_(oji/B,u 2 /B) 2 . Using the 
eigenvectors of A as basis in M. 2 , we obtain a direct sum of two one-dimensional har- 
monic oscillators with frequencies Ba + and Ba_, respectively. This leads to the stated 
form of the eigenvalues. □ 

According to Lemma [!(]] and a simple computation, (JH]) with p 2 = 1 is equivalent to 
J> - 5a + (f, f)(2A; + 1) - Ba_(f, f)(2/ + 1)) + < ^ J> - (2m + l)B)\ 

k,l>0 1 2 m>0 
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with a± given by ( TT5l) . Setting A = fi/B, o~j = Uj/B and a± = a±(a) and substituting 
( fTBT) we can rewrite the desired inequality as 

(A - a+ (2k + 1) - o_ (2/ + 1)) + < (A - (2m + 1)) 2 + , (17) 

k,l>0 + m>0 

and this is what we shall prove. 

3.2. Two inequalities for convex functions. For the proof of (fT7|) we shall need 



Lemma 12. Let <fi be a non-negative convex function on (0, oo) such that J °° (j>(t)dt 
exists. Then 



- - /"OO 

X>(*+±)< / 0(t)dt. 



Proof. Indeed, by the mean- value property of convex functions <p(k + 1) < cfc 
for each k. Now sum over k. □ 

Remark 13. The proof also shows that Yl?k=o + |) — Jo ^(^) ^ ^ or eacn i n teger 
i^. This observation will be useful later. 

The inequality from Lemma [T2l is sufficient to prove a sharp Lieb-Thirring inequality 
in the non-magnetic case, but for the proof of our Theorem[2]we need a more subtle fact 
about convex functions. We note that by the previous lemma hJ^'kLo^fiik + !)) — 
f °° 4>(t) alt for any h > 0. Moreover, ^fclo <HM^ + |)) — ^ Jo^ ^ as ^ ~ ^ by the 
definition of the Riemann integral. The key for proving our sharp result is that, for 
a certain subclass of convex functions, this limit is approached monotonically. More 
precisely, one has 

Lemma 14. Let <fi be a non-negative convex function on (0, oo) such that J °° (p(t) alt 
exists. Assume that <p is differentiate and that (ft 1 is concave. Then the sum 

oo 

& X>(M* + i)) 

k=0 

is decreasing in the parameter h > 0. 
We emphasize that without assumptions on 0' the inequality 

oo oo 

$>(*+!) <hj2 mk+D) 

k=0 k=0 

is not true for all h < 1. Indeed, take for instance <j)(t) = (1 — 1) + and h > 2/3. 

In the proof of this lemma we shall make use of the following well-known fact 
about convex functions: If ^ is a non-negative convex function on (0, oo) such that 
J °° ip(t)dt exists, then ip(t) = / °°(T — t) + dfi(T) for some non-negative measure /i. 
Indeed, it is known that the left-sided derivative d_ip exists everywhere on (0, oo) 
and satisfies ip(b) — ijj(a) = f d^ip{t) dt for < a < b < oo. Moreover, d_ip is 
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increasing and left-continuous, and therefore there is a non-negative measure \x such 
that d-ip(b) — d-ip{a) = //([a, &)). Since lim^oo if)(t) = lim^oo d-ijj(i) = 0, we have 
by Fubini's theorem 



m = - 

as claimed. 



/oo / r \ poo 

I J X[a,oo)(T) dfi{T)J da = J (T — t) + dfi(T) 



Proof. By the fact recalled above (with ip = —(f)') we have 4>{t) = f °°(T — t) 2 + d[i{T) 
for a non-negative measure fi. Hence it suffices to prove the lemma for <p(t) = (T — t) 2 + 
with T > 0. We have to prove that £~ (</>(h(k + \)) + h(k + \)<j)'{h{k + \)) < 0, 
which for our reads 

oo 

{(S -2k- 1) 2 + - 2{2k +\){S-2k- 1)+) < 

k=0 

with S = 2T/h. Choose K G N such that + 1 < S < 2K + 3. Then the left side 
above equals 

K K 

{(S -2k- l) 2 - 2(2k + 1)(S - 2k - 1)) = i s2 ~ AS ( 2k + + 3 ( 2k + X ) 2 ) 

1=0 1=0 

= (K + 1) (S 2 - AS(K + 1) + (2K + 1)(2K + 3)) 
= (K + 1)(S -2K- 1)(S -2K-3). 

This is clearly non-positive for 2K + 1 < S < 2K + 3, thus proving the claim. □ 
3.3. Proof of Theorem [2], We have to prove (fTTj) . By Lemma [T2l for any k 

POO 

(A - a + (2k + 1) - a_(2Z + 1)) + < / (A - a+(2A; + 1) - 2a_t) + 



1 • ^2 



(A - a + (2k + 1)) , . 



4a_ 

A simple computation shows that a + = a + (a) > 1, and hence by Lemma 
a + (A - a + (2k + 1)) 2 + < ^ (A - (2k + l)) 2 , . 

fc>0 fc>0 

The previous two inequalities imply the desired ( TT7|) . □ 

3.4. Proof of Proposition [31 Given < 7 < 1, we want to find u>i = u>2 and -B 

such that the reverse inequality ([7]) holds. We may assume 7 > in the following. 
(The case 7 = can be treated similarly, or one may use the argument of Aizenman 
and Lieb mentioned in the introduction to conclude that a counterexample for 7 = 70 
implies one for all 7 < 70.) 
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By the same computation that lead to ( jTTl) we see that (J7J) can be written as 
(A - a + (2k + 1) - o_ (21 + > 1 fl - £ (A - (2m + 

k,l>0 \l' ) + m > 

with A = fi/B, Oj = Uj/B and a± = a±(a). We will let ui = u 2 and use the notation 
t = a 2 . One can show that a + = 1 + 1 + 0(t 2 ) and a_ = t + 0(t 2 ) as t -> + . We 
now choose A = 3 and recall that a + = a + (a) > 1. This gives us the inequality 

2( 7 + l)a_a + ^ (3 - a+ - a_ (2/ + 1)) \ - 2 7+1 > , 
z>o 

which may be written as 

( 7 + l)al +1 a + J2( x ' l )l ~ 1 > 

z>o 

with x = (3 — a + — a_)/(2a_). Since x = + 0(t)) as t — )■ 0+, we may choose a 
so that x is an integer. In this case we may use the concavity of %p and Remark [TBI to 
bound 

* pX+l/2 -, 

This shows that 

(7 + l)al +1 a + J2( x - l )l > a + {( a - x + «-/2) 7+1 - (a-/2) 7+1 ) 
z>o 

= a + (((3 - a + )/2y+ l - (a_/2)^ +1 ) 

= (1 + 1 + 0(t 2 )) (1 - t/2 + 0(t 2 )Y +1 + 0(^ +1 ) 

= 1 + — p4 + 0(t 7+1 ) . 
Since this is strictly larger than 1 for sufficiently small t, we have proved our claim. □ 
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